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Several years ago it was conjectured in the so-called Roma Approach [Q], that gauge fixing is an 
essential ingredient in the lattice formulation of chiral gauge theories. In this paper we discuss in 
detail how the gauge-fixing approach may be realized. As in the usual (gauge invariant) lattice 
formulation, the continuum limit corresponds to a gaussian fixed point, that now controls both the 
transversal and the longitudinal modes of the gauge field. A key role is played by a new phase 
transition separating a conventional Higgs or Higgs-confinement phase, from a phase with broken 
rotational invariance. In the continuum limit we expect to find a scaling region, where the lattice 
correlators reproduce the euclidean correlation functions of the target (chiral) gauge theory, in the 
corresponding continuum gauge. 
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I. INTRODUCTION 

The great difficulty in constructing chiral gauge theo- 
ries, such as the Standard Model, using the lattice regu- 
larization is related to the doubling problem |^-|^. In lat- 
tice QCD, species doubling occurs when the discretized 
fermion action has an unwanted symmetry that should 
be anomalous in the continuum. When one uses Wil- 
son fermions, the Wilson term eliminates the doublers, 
at the price of breaking all the axial symmetries explic- 
itly. In the continuum limit, one expects to recover the 
axial symmetries, except the anomalous U(l), by tuning 
the fermion hopping parameter to a critical value. 

In the lattice discretization of a chiral gauge theory, 
one has to account for the fact that a Weyl fermion 
in a complex representation contributes to the gauge 
anomaly. This means that a lattice action for a single 
chiral fermion cannot be gauge invariant. We will as- 
sume below that, as in the continuum, the lattice fermion 
action involves a sum over different complex representa- 
tions, whose total gauge anomaly is zero. The question is 
to what extent the violations of gauge invariance, coming 
from the individual representations, cancel each other. 

Consider the regularized effective action obtained by 
integrating out an anomaly-free set of chiral fermions. In 
the continuum, one can use Dimensional Regularization 
to define the effective action for smooth gauge fields that 
vanish rapidly at infinity. Using the freedom to add local 
counter- terms, the violations of gauge invariance are pro- 
portional to the dimensionless parameter e = d — 4, and 
so they vanish in the limit e — > 0. This extends to topo- 
logically non-trivial background fields, using for example 
the (^-function regularization 

On the lattice one encounters a fundamentally differ- 
ent situation. The lattice spacing is a dimensionful pa- 
rameter that plays a dual role. First, it provide a UV 
cutoff, by replacing the infinite range of momentum inte- 
grals with an integration over the periodic Brillouin zone. 
In addition, the lattice spacing enters the (multi-valued) 
mapping from the compact link variables Ux.^, to the Lie- 



algebra valued field. These differences in the global 
structure imply that a generic lattice gauge transforma- 
tion, considered as a mapping that acts on the Fourier 
space of the lattice A^-Held, is qualitatively different from 
the corresponding mapping defined by a continuum gauge 
transformation. The result is that the lattice effective 
action suffers from generic violations of gauge invariance 
which are not controlled by any small parameter. We 
know of no method that ensures the smallness of these 
violations on the entire lattice gauge orbit, at the price 
of tuning any finite number of parameters (see Sect. ^V^ 
or the review article for more details). 

In the so-called Roma Approach ||l|, it was conjectured 
that gauge fixing is a crucial ingredient in the lattice for- 
mulation of chiral gauge theories. A gauge fixing action 
should assign a bigger Boltzmann weight to a smooth 
gauge field, relative to a rough field that belongs to the 
same gauge orbit. This should reduce lattice-artefact vi- 
olations of gauge invariance, because the latter are asso- 
ciated with the roughness of the lattice gauge field. 

In spite of this promising picture, the gauge-fixing ap- 
proach has remained elusive. A naive discretization of the 
Lorentz gauge-fixing action leads to a lattice action that 
has a dense set of lattice Gribov copies with no contin- 
uum counterparts. These lattice artefact Gribov copies 
exist even for the classical vacuum. (Remarkably, the 
proliferation of Gribov copies on the lattice resembles 
the fermion doubling problem in a number of ways.) As 
a result, the Boltzmann weight of too many rough lattice 
configurations is not suppressed. 

In this paper we construct a lattice gauge-fixing ac- 
tion that accommodates this problem. (See ref. |^ for a 
preliminary version of this work.) The gauge- fixing ac- 
tion is associated with a new generic phase diagram. We 
argue that, in this phase diagram, there is a gaussian 
critical point that belongs to the universality class of a 
gauge-fixed continuum theory. In comparison with the 
gauge- invariant lattice definition of QCD, the weak cou- 
pling limit here controls not only the transversal modes 
of the gauge field, but also the longitudinal ones. 
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As discussed above, gauge invariance of the target con- 
tinuum theory, as well as the residual BRST invariance, 
are both explicitly broken on the lattice. By tuning a fi- 
nite number of counter-terms, one hopes to recover BRST 
invariance in the continuum limit, provided the fermion 
spectrum is anomaly-free j^. The BRST identity that 
requires the vanishing of the renormalized gauge boson 
mass, = 0, plays a key role. Since the regularization 
is not gauge invariant, a mass counter-term has to be in- 
troduced, and its parameter needs to be tuned, in order 
to enforce this BRST identity. Usually, a negative renor- 
malized mass-squared indicates spontaneous symmetry 
breaking. Here we encounter a new feature, namely, the 
gauge field condenses if its mass-squared parameter be- 
comes too negative. This implies that the new critical 
point is located on the boundary between a conventional 
phase, which is invariant under lattice rotations, and a 
new phase where the lattice rotation symmetry is broken 
spontaneously by the VEV of the gauge field. 

The construction of the gauge-fixing action is pre- 
sented in Sect. ||. The main results are (a) the gauge- 
fixing action has a unique absolute minimum, Ux^fj, = I, 
and (b) perturbation theory around this minimum is 
manifestly renormalizable. In Sect. Ill we discuss a sim- 



ple chiral fermion action. The validity of perturbation 
theory implies the onset of a scaling behaviour in the 
weak-coupling limit. Up to the regularization-dcpcndcnt 
counter-terms, the continuum lagrangian that controls 
the scaling behaviour can be read off from the marginal 
and relevant terms of the lattice action. The scaling re- 
gion should therefore faithfully reproduce the correlation 
functions of the target chiral gauge theory, in the corre- 
sponding gauge. 

Without the new gauge-fixing action, the fermion ac- 
tion of Sect. Ill does not lead to a chiral gauge theory in 
the continuum limit. The longitudinal modes fluctuate 
strongly, and their non-perturbative dynamics ultimately 
renders the fermion spectrum vector-like. (This applies 
to many other chiral fermion proposals, see ref. [[7|j9[|.) In 
the second part of this paper, we examine the dynamics 
of the lattice longitudinal modes from a broader point of 
view. We explain the problems created by this dynamics, 
and how they may be solved within the present approach. 

In Sect. IV we discuss the lattice effective action, and 
the role of lattice artefact Gribov copies. In Sect. we 
discuss the complete phase diagram in the limit of a van- 
ishing gauge coupling. In Sect. VI we explain how our 



approach evades the No-Go theorems. Several open ques- 
tions arc discus sed i n Sect. |VII|, and our conclusions are 
offered in Sect. 



VIII 



II. CONSTRUCTING THE LATTICE 
GAUGE-FIXING ACTION 

A. The phase transition associated with 
a critical vector boson 

This section is devoted to a step-by-step construction 
of the gauge-fixing action. As discussed in the introduc- 
tion, the lattice-regularized theory has no symmetry that 
protects the masslessness of the vector bosons. There- 
fore, in the relevant part of the phase diagram, the lattice 
vector field is generically not critical. Now, according to 
the standard lore, the correlation length should diverge 
close to a continuous phase transition associated with the 
condensation a Bose field. In this paper we wish to apply 
this to the lattice vector field. 

As a preliminary requirement for a continuous tran- 
sition, one needs a higher-power term, that stabilizes 
the classical potential Vc\{Afj_) when the coefficient of the 
quadratic term changes sign. In a weakly-coupled theory, 
the actual location of the transition should be close to its 
tree-level value, and near the transition one expects the 
onset of a scaling behaviour governed by renormalized 
(continuum) perturbation theory. 

Our aim is to achieve criticality of the lattice vector 
field, very much like the way this is done in the famil- 
iar $^ theory. However, going from a spin-0 to a spin-1 
field presents new difficulties. The lattice theory is for- 
mulated in terms of the link variables Ux,^,, which are 
group-valued parallel transporters. On the other hand, 
renormalized perturbation theory, that governs the scal- 
ing region, is more naturally formulated in terms of the 
Lie-algebra valued field. Thus, it takes some trial and 
error to find the [/^-dependent action that best suits our 
purpose. 

Another complication arises because not every renor- 
malizable vector theory is unitary. A unitary, physical 
Hilbert space exists if and only if the vector theory is 
actually an anomaly-free gauge theory in a gauge fixed 
form. This requires us to choose the lattice action, such 
that the marginal gauge symmetry breaking terms in the 
tree-level vector lagrangian have the form 
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(gauge condition)^ 



(2.1) 



An appropriate Faddeev-Popov ghost action will be nec- 
essary too. 



B. A higher-derivative Higgs action 

Our starting point is the lattice action 

S^Sg{U) + Sh{<I>,U). (2.2) 

Here Sg{U) is the usual plaquette action. The Higgs 
action is 
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5H = tr^(-K0tn((7)0 + /i0^n2(t/)0) , (2.3) 



(2.6) 



where 



(2.4) 



is the standard nearest-neighbour covariant laplacian. 
The lattice spacing a is equal to one. Both Ux,^i and (jjx 
take v alue s in a Lie group G. The first term on the r.h.s. 
of eq. (2^) is a conventional lattice Higgs action, whereas 
the second term is a higher derivative (HD) action. 

(HD actions were recently discussed by Jansen, Kuti 
and Liu . Here we are interested in a different critical 
point from the one studied in ref. At the techni- 

cal level, this allows us to introduced only a laplacian- 
squared HD term, whereas for the purpose of ref. it 
was crucial to introduce also a laplacian-cubed one.) 

The action eq. ( |2.2| ) is gauge invariant, where the lat- 
tice gauge transformation is given by Ux,ti 9xUx,fj.gl+j^ 
and ipx gx4>x for gx S G. Now, since 4>x & G too, we 
may use the lattice gauge invariance to completely elim- 
inate the (px field. Note that this operation affects only 
Sh- We introduce the notation 



Sv{U)^SH{(b, U) 



(2.5) 



The subscript of Sy stands for "vector". Sh can be 
recovered from Sy by making the substitution Ux,fi — > 
4>l:Ux,fj.4>x+p,- (The significance of the (px field, which is 
associated with the longitudinal degrees of freedom, is 
discussed in Sect. ^ and Sect. f7\ ) 

W e w ill denote the first formulation of the theory 
(eq. ( |2.2| )) as the Higgs picture. The alternative formula- 
tion (eq. ( |2.5| )) where only Ux^fj, (but not (j>x) is present, 
is called the vector picture. The equality of the partition 
functions in the two pictures extends to observables. Any 
observable in the vector picture is mapped to a gauge in- 
variant observable in Higgs picture, and vice versa. Thus, 
we are dealing with two mathematically equivalent for- 
mulations of the same theory jll]] . 

In this section we assume k ^ 1. The physics in this 
parameter range is more easily accounted for in the vec- 
tor picture, which is used below to study the classical 
potential, and to set up the weak-coupling expansion. In 
the vector picture, the gauge non-invariance of the action 
resides in Sy- Therefore, we will ultimately demand that 
the marginal terms in Sy have the form of a gauge fixing 
action, cf. (pT]). 



and considering a constant Af^ field, the action Sy leads 
to the following classical potential 



Vei = nTigoA^) + kT^igoA^) 



J'igoA^) = 2^ (l - cos(.goA^) 



(2.7) 
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Note that Sg is zero for a constant abelian field. 

For K > 0, the absolute minimum of the classical po- 
tential is A^ = (mod 27r/go)- Since the quadratic term 
in Vci comes only from the K-term, a non-zero vector 
condensate arises for k < 0. The classical features of the 
transition can be determined by keeping only the lead- 
ing term in the expansion of J-'{goAfj,), separately for the 
K- and K-terms. (This approximation is consistent for 
|k| ^ K.) The result is the quartic potential 
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(2.9) 



Eq. (2^) closely resembles the potential of a theory. 
For small negative k, the minimum is 



1 

go 



K<0. 



(2.10) 



Eq. ( 2.10| ) exhibits the mean-field critical exponent 1/2. 
It is easy to check that this is the absolute minimum of 
the classical potential. 



We note that eq. (2.10) is invariant under arbitrary 
S0(4) rotations, refiecting the symmetry of the approxi- 
mate potential eq. ( |2.9| ). When higher-order corrections 
are taken into account, the rotational symmetry of the 
potential is reduced to the lattice hypercubic symmetry. 

Below, the phase with a non-zero vector condensate 
will be denoted as the FMD phase. We will speak about 
the FMD transition, referring to the transition from the 
rotationally invariant phase to the FMD phase in the 
large-K region. FMD stands for ferromagnetic direc- 
tional. The preferred spacetime direction of the FMD 
phase is defined by the vectorial VEV. For gQ ^ 0, there 
are no Goldstone bosons in the FMD phase, because the 
lattice rotation group is discrete. The limiting go = the- 
ory is discussed in Sect. and in particular we explain 
there in what sense the FMD phase is ferromagnetic. 



C. The classical potential 

For simplicity we consider the classical potential in the 
U(l) case. The essential features generalize to the non- 
abelian case. Making use of the standard weak-coupling 
expansion 



D. The weak-coupling expansion 

We now want to study fluctuations around the clas- 
sical vacuum Vx^p. — I (equivalently A^ = 0) in the 
rotationally-invariant phase, close to the FMD transition 
where the theory defined by eqs. (2-5) is expected to be 
critical. The FMD transition is given by k = in the 
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classical approximation. As mentioned earlier, we are 
assuming k ^ 1. We therefore focus on the HD term 
in its vector picture form. Relaxing the assumption of a 
constant field, we find 
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\ A' 
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(2.11) 



where the dots stand for irrelevant operators. is the 
backward lattice derivative, defined as A^/j; = fx^ fx-ji 
for any func tion fx- 

Eq. (2.11) contains a longitudinal kinetic term. We 
define 
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= f^9o : 



(2.12) 



and we will assume that is an 0(1) parameter. This 
means that the longitudinal kinetic term belongs to the 
tree-level lagrangian. Remember that a transversal ki- 
netic term is provided by the gauge invariant plaquette 
action Sq- Finally, we assume that the tree- level vec- 
tor boson mass is zero. Under these assumptions, the 
tree-level vector propagator is 



(2.13) 



marginal term - the longitudinal kinetic term - has the 
form of a gauge-fixing action for the linear Lorentz gauge 
d ■ A — 0. Alternatively, the new HD term can lead to a 
mixed marginal term proportional to {d ■ A)A'^. In this 
case one recovers the non-linear gauge d ■ A + gA^ = 0. 

The linear gauge d ■ A — is more familiar, and less 
complicated to implement in perturbation theory. More- 
over, the above non-linear gauge is consistent only for 
U(l) or SU(N)xU(l), whereas the linear gauge is con- 
sistent for any gauge group. The linear gauge has, how- 
ever, one technical disadvantage. The quartic term in 
eq. (2.11) is the stabilizing term of the classical poten- 
tial (see eq. (pj|)). In its absence, one has to reanalyze 
the classical potential, and make sure that it is stabilized 
by a higher-power term (in practice this is an A^ term). 
This task is done in ref. , which is henceforth referred 
to as //. 

Here we will consider only the non-linear gauge. Since 
the necessary mixed term contains a derivative, one can 
modify the HD action while leaving the classical potential 
intact. This simplifies our task, as the large-K study of 
the phase diagram in Sect. II C remains valid. The new 
HD action is 



Qll.l. 
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■ tr 



(^0^n'(c/)</> + 2B^A;v;j 



(2.16) 



where 



(2.14) 
(2.15) 



and = 2sin(pp/2). 

Massless weak-coupling per turba tion theory is defined 
by the vector propagator eq. ( 2.13 ), and by a set of ver- 
tices which can be read off from the lattice action using 
eq. (O) in the usual way. 



E. The gauge-fixing action 

In view of the presence of kinetic terms for all polariza- 
tions, lattice perturbation theory is manifestly renormal- 
izable. According to the standard lore, renormalizabil- 
ity implies a Lorentz invariant scaling behaviour in the 
vicinity of the gaussian critical point go = 1/k = 0. The 
scaling behaviour is achieved by tuning a finite number 
of counter-terms, that correspond to the relevant and the 
marginal operators. 

At this stage, the marginal gauge symm etry b reaking 
terms in the tree- level vector action (see eq. ( 2.11 )) do not 
have the form of a gauge- fixing action, cf. ( 2.1 ). The way 
to remedy this is to add another term to the HD action. 
There are two options. The new term can be chosen 



to cancel the quartic term in eq. (2.11). The remaining 



Vx. 



E 



Vx- 



V, 



X,fJ. 



2i 



7 {4>\Ux,^,4>x+^i - h-c.) 



(2.17) 



(2.18) 



Going to the vector picture and applying the weak- 
coupling expansion, we have 



(50^, 



(2.19) 



It is easy to check that the desired mixed term is now 
present in the tree-level vector action. We comment that, 
in the Higgs picture, is a gauge-invariant local vector 
field, whose expectation value serves as an order param- 
eter for the FMD phase. (The corresponding order pa- 
rameter in the vector picture is the expectation value of 

(V^> = 50 {A^), where 



=/■) classical level, 



the latter is given by eq. ( ^.10 ). 

For the laplacian-squared HD action (see eq. ( ^.3[ )), 
it is evident that Ux,^, = I is the unique absolute min- 
imum for all configurations, and not only for the con- 
stant ones considered in the classical potential (we as- 
sume (px — I)- This property, which is necessary to val- 
idate the weak-coupling expansion, applies to the new 
HD action ( 2.16| ) as well. The proof is given in //. 
The symmetric combination used in the definition of Bx 
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(eq. ( ^.17 )), which does not affect the marginal term con- 
tained in BJ2^ essential for the proof. 
We define the lattice gauge-fixing action to be 



As expected, S'^^^' has the classical continuum limit 
(l/2^o)(5 • A + gJ^Y . Because of the irrelevant terms it 
contains, one cannot write S'^^' as the (sum over x of the) 
square of a local function of the L/^-s. Consequently, the 
gauge-fixed lattice action is not invariant under BRST 
transformations. 

It is interesting that the breaking of (gauge and) BRST 
invariance is a common feature of the chiral fermion ac- 
tion and the gauge- fixing action. In the case of the gauge- 
fixing action, it has to be so because of a theorem by 
Neuberger |Q, which asserts that any lattice BRST- 
invariant (gauge-fixed) partition function must vanish 
due to lattice artefact Gribov copies. We return to the 
role of lattice Gribov copies in Sect. 

Before we introduce fermions, the complete lattice ac- 
tion (in the vector picture) is therefore 



<5g + 



n.l. 



5n.l. 
fp 



-Sc. 



(2.21) 



For the non-linear gauge, the continuum Faddeev-Popov 
action involves the operator (we suppress the group struc- 
ture constants) 9^ -f igA ■ d + g{A, d + igA}. The last 
term is absent in the case of the linear gauge. For the 
discretization of we take the standard (free) lattice 
laplacian. We are discretizing a second order operator, 
and our choice avoids the appearance of any FP dou- 
blers. For the interaction terms, any lattice operator 
with the correct classical continuum limit should do. For 
the discretization of gA ■ d, for example, one can take 
Vp^/S.^^ where is the antisymmetric difference op- 
erator. (Since BRST symmetry is broken anyway by the 
gauge-fixing action, wc make no attempt to preserve any 
exact relation between the discretized versions of and 
A ■ d.) We note that the ghost fields contribute to the 
effective potential only through loops, and so they do not 
modify the tree-level considerations. 

Set is the counter-term action. The role of Set is to 
enforce BRST invariance in the low momentum limit of 
lattice perturbation theory The BRST symmetry 
is violated in particular by (marginal) S0(4)-breaking 
lattice operators. Therefore, enforcing BRST invariance 
should also restore full S0(4)-invariance in the contin- 
uum limit. The counter-term action is more naturally 
written in terms of A^. We define Set as a local func- 



tional of the J7^-s by trading A^ with using eq. ( 2.19 ). 
(The second term in the expansion of V^, which breaks 
SO (4) invariance, is needed only for the dimension-two 
mass counter-term. In all other cases one simply replaces 
goAf^ with V^.) 

The BRST identity = 0, which says that the (renor- 
malized) vector boson mass must vanish to all orders in 



perturbation theory, is consistent with taking the contin- 
uum limit at the FMD transition. As a mass counter- 
term one can take the K-term in eq. ( ^.3| ). This means 
that K, is tuned to Kd.idOj Co)i where in perturbation the- 
ory Kc.i. = J2n>i Cn(^o)ffo^"~^''- Notc that the coefficient 
of the mass term in eq. (2.E) is Kg^. The absence of an 



0(l/gg) term in the expansion of Kd. is in agreement 
with the vanishing of the tree-level vector boson mass. 

In this paper we have simplified things by consid- 
ering only the most important counter-term, namely, 
the dimension-two mass term. In the case of the non- 
linear gauge, the next most important counter-term is 
the dimension-four S0(4)-breaking term J^fj.^t-' 
the linear gauge, the classical potential is stabilized by 
an A^ term, and (^ A'j^)'^ too occurs only as a counter- 
term. In this case, the effect of the dimension-four non- 
derivative counter-terms is discussed in //. One finds 
that the conventional (hyper-cubic invariant) phase and 
the FMD phase both extend into the higher-dimensional 
phase diagram. Also, the FMD transition remains con- 
tinuous when the dimension- four counter-terms are tuned 
to their critical values. The crucial features leading to 
these conclusions are (a) it is justified to expand up 
to a finite order in (equal to the dimension of the sta- 
bilizing term) when looking for the absolute minimum of 
the potential; (b) the coefficients of the counter-terms are 
0(1), whereas the coefficient of the gauge-fixing action is 
O(l/go)- Since these features are true in the case of the 
non-linear gauge as well, we expect a similar robustness 
against the inclusion of additional counter-terms. 

In this section we have discussed the phase diagram 
only in the large-K limit. The phas e diagram for arbi- 
trary K and K is studied in Sect. |V A . This study , as well 
as additional arguments presented in Sect. VIA , further 
clarify why the continuum limit of the gauge-fixing ap- 
proach should be defined at the FMD transition. 



III. CHIRAL FERMIONS 

In a gauge invariant lattice theory, the minimum of the 
plaquette action is unique up to a gauge transformation, 
and the transversal kinetic term is sufficient to define a 
valid weak-coupling expansion. In the absence of gauge 
invariance, there exists a valid weak-coupl ing e xpansion 
provided the gauge- fixing action of Sect. HE is added 
to the plaquette action. This applies also to the gauge- 
fixing action presented in II for the linear gauge d-A = 0. 
Using either of these gauge-fixing actions, there is a lot 
of freedom in the choice of the chiral fermion action. We 
consider here (in the vector picture) an action which is 
the most economic in the number of fermionic degrees of 
freedom (For related work see ref. |l^Jl^.) Other 

fermion actions have certain advantages over the one pre- 
sented here, and in particular they can reduce the re- 
quired fine-tuning. 

According to ref. [ p^ , one introduces a two-component 
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lattice fermion field Xx^ to account for a single Weyl 
fermion in the target continuum theory. The fermion 
action is (suppressing coordinates summations) 



SF^y,X^^.D^,{U)x-J{xUx + '^■C■), (3.1) 



Dxy,,,{U) ^ ^(4+A,J/C^x.p - Sx-il,yUl^) . (3.2) 



Here d^y is the free lattice laplacian (eq. (2^) for Ux,fj, — 

I), and XxXy = ^a0Xx,aXy,f3 = xI<^Xy where e is the an- 
tisymmetric two-by-t wo m atrix. We assume w = 0(1). 
The first term in eq. ( ^.l| ) is the naive lattice discretiza- 
tion of the continuum Weyl action. The second term is 
a Majorana- Wilson (MW) term, that breaks explicitly 
gauge invariance as well as the fermion number symme- 
try. The latter is unwanted, because fermion number is 
not conserved in the continuum theory. 

In order to understand the properties of the lattice 
fermion path integr al, it is convenient to recast the 
fermion action eq. (3.1) in terms of four components 
il)M- By definition, Pl^Pm — X and 
where Pr^l = 3(1 ± 75) denote chiral- 
ipM is not an independent field, and is 



and 

T 



fields ■0M 
Pr^^m = ex 
ity projectors 
given by 



M 



(3.3) 



Here C is the antisymmetric four-by-four charge conju- 
gation matrix, obeying — —1, 7JC ~ —C^^^, and 
75C = C75. In terms of these four component fields, the 
fermion action takes the form 



= 9 E \i^D^{U)Pl + ^^D^{U*)Pe 



M 



^ 7 

-—WmUWm ■ 



(3.4) 



Let us first examine eq. (3.4) in perturbation theory. 
The tree-level fermion propagator is the (massless) Wil- 
son propagator p4| . Because of the Majorana-like condi- 
tion eq. (3.3), the symmetry factors in Feynman graphs 
are the same as for Majorana fermions. Now, if we go to 
the small momentum limit, we find that the two chiral- 
ities of ipM couple to the gauge field according to com- 
plex conjugate representati ons of the gauge group. One 
sees that the role of eq. (3^) is to consistently main- 
tain the identification Pl^^m Xi PbMim Xj ^it the 
level of Feynman diagrams. That eq. (3^) correctly de- 
scribes a single left-handed Weyl fermion, can be veri- 
fied by calculating the non-analytic part of one-fermion- 
loop diagrams, that should agree with the continuum 
result in the limit of a vanishing external momentum 
(the role of counter-terms is discussed below). As usual, 
the non-analytic contribution comes from an infinitesi- 
mal neighbourhood of the origin in the Brillouin zone. 
In this neighbourhood one can neglect the Wilson term 



in both numerators and denominators. The left-handed 
and right-handed components of ipM are no longer cou- 
pled, and one can reexpress the Feynman integrand in 
terms of the continuum propagator for a single Weyl 
fermion. (A left-handed fermion loop is equal to a right- 
handed fermion loop in the complex conjugate represen- 
tation. The "double-counting" is compensated by a one- 
half symmetry factor for each closed fermion loop, which 
arises from the Majorana-like condition eq. ( |3.3D .) 

We next discuss the rigorous definition of the fermionic 
path integral. We introduce the 2N x 27V fermion ma- 
trix Qaf3 by writing the action in the generic form Sp — 
\ 13 V'm Qap V'm- (The chargc conjugation matrix C 
is absorbed into the definition of Q.) It is easy to check 
that Q is antisymmetric. The fermion path integral takes 
the following form [O 



/ n exp j ^ XI V' 



= pf 



(3.5) 



There is no integration over the (dependent) variables 
•ipM- According to eq. ( ^.5[ ), the fermionic path inte- 
gral is a Pfaffian. In general, pf((5/2) is complex, as 
expected from the euclidean path integral for a single 
Weyl fermion. 

(As a further check that our fermion path integral 
describes a Weyl fermion, we can consider a "two- 
generation" model, where each complex representation 
occurs twice in the fermion spectrum. Using the iden- 
tity pf^((5/2) = det(Q), that holds for a general anti- 
symmeric matrix, this two-generati on m odel can be de- 
fined by an action similar to eq. (3.4), where we now 
drop an overall one- half factor, substitute Va/ ~^ V'Dj 
V'Af — > V'Dj and regard and ipD as independent Dirac- 
like variables. The counting of degrees of freedom is now 
straightforward. Since the two chiralities oi ipo belong 
to complex conjugate representations, this action actu- 
ally describes two left-handed Weyl fermions in the same 
complex representation.) 

Until now we have implicitly discussed the fermions 
in the background of a fixed external gauge field. The 
main result of the previous section is that, with the 
gauge-fixing action, perturbation theory is valid for a dy- 
namical gauge field as well. Therefore, with appropriate 
counter-terms, the continuum fields describing the scal- 
ing behaviour are in one-to-one correspondence with the 
massless poles of the various tree-level propagators. If we 
choose an anomaly-free fermion spectrum, the scaling re- 
gion should be governed by a continuum chiral gauge the- 
ory, in the relevant gauge. We note that if one chooses an 
anomalous fermion spectrum, the scaling region will still 
be governed by a renormalizable lagrangian, but BRST 
invariance and, hence, unitarity will be violated. 
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Finally, let us discuss the fermion mass counter-terms. 
As with ordinary Wilson fermions, a mass counter-term 
mo^Mi'M is necessary to maintain the masslessness of 
each chiral fermion. (A different fermion action that 
does not require mass counter-terms will be discussed 
elsewhere |36 .) The renormalized Majorana-like mass 
is proportional to (mg — rric), where rric is (minus) the 
fluctuations- induced mass. If (mo — TOc) is small but non- 
zero, BRST invariance will be explicitly broken in the 
scaling region. The scaling behaviour is then governed 
by a renormalizable continuum theory which is not gauge 
invariant (hence also non- unitary). By tuning toq to mc, 
assuming all other counter-terms already have their crit- 
ical values, we recover BRST invariance simultaneously 
with the masslessness of the chiral fermions. (The situa- 
tion on the lattice is similar to what one would encounter 
in the continuum, if a gauge non-invariant regularization 
is employed for a chiral gauge theory. As on the lat- 
tice, Majorana-like mass counter-terms may be needed, 
alongside with other gauge non-invariant counter-terms, 
to cancel the breaking of gauge invariance induced by 
the regularization, and to ensure that the renormalized 
amplitudes are gauge invariant.) 

An important question in the literature on lattice chi- 
ral gauge theories, is how to correctly reproduce fermion 
number violation in the continuum limit. Different so- 
lutions have been proposed to the problem |2l|. We 
hope that the present approach can shed new light on it. 

For definiteness, we adopt the strategy of rcf. p4[ . 
Namely, we demand that the lattice fermion action 
should have no symmetry which is not present in t he t ar- 
get continuum theory. Now, while the action eq. ( |3.l| ) is 
not invariant under global U(l) transformations with an 
arbitrary phase, it is still invariant under the residual dis- 
crete symmetry x ~^ ~Xi X ~^ ~X- This symmetry im- 
plies a (mod 2) conservation law for each fermion species, 
which still causes a problem. Consider for definiteness an 
SU(5) GUT, with one generation that contains a 5 and 
a 10. In an instanton background, the numbers of zero 
modes for these representations are respectively one and 
three. This is in conflict with the above (mod 2) con- 
servation laws. Thus, on top of the MW terms present 



in eq. (3.1) for each representation, one has to introduce 
an additional gauge-noninvariant MW term that couples 
the 5 and the 10. (As a result, a Majorana-like mass 
counter-term that mixes the 5 and the 10 will be nec- 
essary too.) With this new MW term, the remaining 
discrete symmetry leads only to (mod 2) conservation of 
the total fermion number for each generation. 



IV. WHY GAUGE FIXING 

We now return to the lattice effective action, and con- 
sider some of its properties in more detail. We keep the 
discussion at an informal level. Our approach has been 
presented in detail in the previous sections, and the aim 



here is to clarify the nature of the problems that it is 
meant to solve. 

In this section we assume that the lattice chiral fermion 
action Sp is bilinear in the fermion fields, and that it 
depends in addition only on the link variables Ux,^,- (This 
corresponds to the vector picture.) It is also assumed 
that Sp is (mildly) local. The lattice spacing a will be 
shown explicitly in this section. 

The difficulties encountered in the construction of lat- 
tice chiral gauge theories can be addressed at a more rig- 
orous level 1^,0. This complementary discussion, which 
focuses on the robustness of the Nielsen-Ninomiya theo- 
rem, is given in Sect. VI. 



A. Rough lattice gauge transformations 
and the need for gauge fixing 

The lattice effective action is defined by integrating out 
the fermions 

S,s{U) = - log / V^Vij e-^^^^'^'^^ . (4.1) 



Clearly, the well-defined object is exp(— S'off), rather than 
ScS itself. For our purpose it will be sufficient to consider 
the perturbative effective action, and so we will ignore the 
problems associated with the global definition of Scft- 

The variation of SeS, in response to an infinitesimal 
lattice gauge transformation at the point xq, has the fol- 
lowing general form 



SxgScS ~ Co Ox°" 



n>l 



(4.2) 



The « sign indicates that the r.h.s. is computed pertur- 
batively. Note that the gauge field is external, and so 
the gauge-field action is not needed at this stage. O™" 
and 0"' ' ar e local lattice operators that depend on A^, 
cf. eq. (2.6). The dimension of O"'* is 4 -I- n, and the i- 



summation is over linearly independent operators of this 
dimension. O™" is some discretized version of the consis- 
tent anomaly. We assume that all operators of dimension 
less than or equal to four, other than 0^°", have been 
cancelled by counter-terms. If, moreover, we choose a 
set of complex representations that satisfies the anomaly 
cancellation condition, then cq = 0. 



The infinite sum on the r.h.s. of eq. (4.2) accounts 
for lattice artifact violations of gauge invariance. The 
precise form of these violations is model dependent, but 
their existence is generic. As can be easily seen by go- 
ing to momentum space, which is the usual setting for a 
perturbative computation, this sum represents a double 
expansion in |goa^/j| and |ap^|. 

Let us first consider a smooth gauge field A° which is 
characterized by some physical scale Aphys ^ a~^, and 
the correspon ding configuration of link variables ^ de- 
fined via eq. ( p.q ). Since the dimensionful quantities A^^ 
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and Pfj, are O(Aphys), b oth of the above expansion pa- 
rameters are smalL Eq. (^^) is the gradient of SeS with 
respect to a motion inside the lattice gauge orbit. There- 
fore, ScS is (approximately) constant on the orbit in the 
vicinity ^. The constancy of SeS extends to that por- 
tion of the orbit which is reachable from ^ by a smooth 
gauge transformation. 

The problem is that, on the lattice, smooth gauge 
transformations represent a tiny part of the local gauge 



group. Let Ux^^l 



be a generic configuration 



in the orbit of U^^^, and let a'^^ be related to C/i^^l via 
eq. (|2.6|) . Since the operators that occur on the r.h.s. of 
eq. ( |4.2D are not gauge-invariant, they are sensitive to the 
value of gx for x in the vicinity of xq. Now, the ^^.-s on 
different sites a uncorrelated. As a result, the expansion 
parameters l^oa^A*'^! and {aA^A^f^ /A^'ff^ \ are 0(1) for a 
generic lattice gauge transformation. We conclude that 
ScS fails to be (approximately) constant on most of the 
lattice gauge orbit. This is true for any orbit, including 
orbits that have a smooth representative. 



B. Proliferation of lattice Gribov copies 

The above problem stems from the roughness of generic 
lattice gauge transformations. Following ref. [0 we make 
no attempt to reduce the violations of gauge invariance 
at the level of the effective action. Instead, our aim is to 
suppress the Boltzmann weight of rough gauge field con- 
figurations (relative to smooth configurations that belong 
to the same orbit) consistently with the gauge invariance 
of the physical Hilbert space, namely, via gauge fixing. 

In lattice QCD, gauge fixing is a matter of choice, since 
it has no effect on the gauge-invariant observables. Here, 
the fermion action is not gauge invariant. As a result, 
the gauge-fixing method is an integral part of the defini- 
tion of the theory. Different gauge-fixing methods may 
in general give rise to different phase diagrams with dif- 
ferent critical points. There is no guarantee that every 
gauge-fixing method will lead to a non-trivial continuum 
limit, let alone to a chiral gauge theory. 

Still, in order to make progress, one has to choose some 
gauge-fixing method. Vink proposed to use the lapla- 
cian gauge, where a maximally-smooth representative is 
chosen on each gauge orbit by global minimization. The 
laplacian gauge is highly non-local, and this creates dif- 
ficulties both in the analytic and in the numerical study 
of this method. 

As discussed in detail in Sect. ||, we build a local 
gauge-fixing lattice action, that (a) has the unique ab- 
solute minimum Ux,fi = I, and (b) reduces to a covari- 
ant gauge-fixing action in the classical continuum limit. 
Our gauge-fixing actions (eqs. (16-18) for the non-linear 
gauge, see // for the linear gauge) are clearly not the 
most naive discretizations of the corresponding contin- 
uum actions. Focusing for simplicity on the linear gauge 



d- A = 0, let us examine what goes wrong with a naively- 
discretized gauge-fixing action. We thus consider the fol- 
lowing action 



cnaivc 



1 



2^030 



tr 



(4.3) 



(4.4) 



Note that g^^ A~Vx^^ reduces in the classical contin- 
uum limit to d ■ A, as it should. 

What is common to our gauge- fix ing action(s) and 
to the naive gauge-fixing action eq. (4.3), is that they 
contain a longitudinal kinetic term. The trouble with 
gnanc ^j^^^ supports a dcnsc set of Gribov copies 
for the identity field Ux,fj. = I- Each of these Gribov 
copies is a classical vacuum of So + -Sgl"™. The su- 
perposition of contributions coming from all these clas- 
sical vacua, which cannot be calculated perturbatively, 
may ultimately render the fermion spectrum vector-like. 
Even without fermions, stability of the classical potential 
is lost, and it is unclear whether the FM-FMD transi- 
tion (associated with a divergent vector-field correlation 
length) can be maintained in the weak-coupling limit. (If 
the gauge-fixing action in eq. (2.21) is replaced by 5""^' ™, 
the resulting classical potential is identically zero. If we 
assume k = 0(l/gg), which implies the presence of a 
tree-level mass term, one has y^."^'™ = K!F{goAf^) (com- 
pare eq. i^)). The minimum of FJJ^'™ is = for 
K > 0, and A^ = it/ go for k < 0.) 

We now demonstrate the existence of a dense set of 
lattice Gribov copies for the identity field. Consider first 
the U(l) case. The condition tJa; = is satisfied if the 
imaginary part of Ux.fi is zero everywhere. The latter 
is true if we consider only lattice gauge transformations 
where g^; = ±1. In other words, in spite of the presence 
of the gauge- fixing action .Sgf ™, the Gribov copies of the 
identity field still exhibit a local Z2 symmetry. (This is 
also true for the non-linear gauge, if one replaces Gx h^ 
eq. dij) by ^-1- = E^(A;:^.,M + Vl^).) An "elemen- 
tary" Gribov copy is created if we choose gx = —I for 
X = Xq, and gx = ^ elsewhere. This clearly shows that 
the Gribov copies are local lattice artefacts. 

The above example generalizes to non-abelian groups. 
In the case of SU(2N) and S0(2N) groups, simply re- 
place ±1 by ±/. Moreover, for any SU(N) group, one 
can choose an SU(2) subgroup (which for simplicity we 
assume to lie at the top left corner) and repeat the above 
construction with gx = diag{±l, ±1,1,1,...). The dis- 
crete local symmetry of the Gribov copies is therefore 
larger than Z2 in the general case. For a more detailed 
discussion of lattice Gribov copies see ref. . 

A number of remarkable similarities draw us to say 
that the proliferation of Gribov copies is the spin-1 
counter-part of the fermion doubling problem. In both 
cases, one deals with the discretization of a first-order 



8 



differential operator: in the spin-i case, this is the Dirac 
(or Weyl) equation; in the spi n-1 c ase, this is a covari- 
ant gauge condition (see Sect. HE). In both cases the 
problem arises when a non-compact continuum variable 
is replaced with a compact lattice variable: fcrmion dou- 
bling arises because, unlike in the continuum, the lattice 
momentum is periodic; in the spin-1 case, also the non- 
compact continuum gauge field is replaced with compact 
group variables. In both cases, there are theorems that 
establish an impasse under certain mild-looking condi- 
tions: the Karsten-Smit and Nielsen-Ninomiya Q the- 
orems which predict fermion doubling, and Neuberger's 
theorem which asserts that any BRST-invariant par- 
tition function must vanish identically. Finally, in both 
cases the solution is to reduce the symmetry of the lat- 
tice theory, by adding irrelevant terms to the naively- 
discretized action. In the case of Wilson fermions, this 
is the role of the Wilson term. As for our gauge-fixing 
action(s), one can show (see //) that it is equal to the 
square of a discretized gauge condition, plus irrelevant 
terms. Thus, again, the irrelevant terms reduce the sym- 
metry, this time by breaking explicitly BRST invariance. 



C. Other approaches 

A different approach to the dynamical problems cre- 
ated by rough lattice gauge transformations is to adopt a 
more sophisticated definition for the effective action. The 
prominent representatives of this approach are the inter- 
polation method |^,^ and the overlap formalism p^ . 

In the interpolation method, one constructs a contin- 
uum interpolating field = ^;^°"*(y; C/^,,^), y € i?^, 
for each configuration of the lattice gauge field. The de- 
terminant of the Weyl operator, a ■ {d + igA'^°'^^), is then 
defined using a separate regulator. Associated with the 
fermion regularization is a new cutoff parameter, denoted 
generically At, which must be sent to infinity before the 
lattice spacing a is sent to zero. (A concrete method ||2^] 
is to discretize the Weyl action on a finer lattice with a 
lattice spacin g Of <C a, using e.g. the fermion action of 
ref (eq. (|^). In this case Af = a^^^) 

Consistent regularizations of the Weyl determinant 
break gauge invariance for finite values of the cutoff, 
even when the fermion spectrum is anomaly-free. There- 
fore, gauge-noninvariant counter-terms are needed in the 
interpolation method too. (It has been proposed that 
gauge-noninvariant counter-terms may be avoided, if the 
real part of the effective action is regulated separately 
from the imaginary part ||2^ . In spite of attempts in this 
direction ]25[ |, it remains unclear whether this procedure 
can be implemented beyond perturbation theory without 
violating locality, and, eventually, unitarity.) 

For given Ux,^-s, the interpolating field assigns a local 
winding number to each hypercube (in the non-abelian 
case), or to each plaquette (in the abelian case). In the 
non-abelian case, this is the winding number of the con- 



tinuum gauge transformation defined on the faces of the 
hypercube, that brings the interpolating field to a pre- 
scribed axial gauge; in the abelian case, the continuum 
gauge transformation is defined on the perimeter of each 
plaquette. Now, a fundamental requirement is that the 
fermion determinant should be gauge invariant in the 
limit At — > cx). Gauge invariance can be established only 
if yl™"*(y) is globally bounded [Q|5). Gauge- invariance 
is therefore recovered in the limit Af ~> oo only on that 
portion of the lattice gauge orbit, where all the local 
winding numbers are zero [Q. The solution is to ap- 
ply a gauge transformation that sets all local winding 
numbers to zero before computing the fermion determi- 
nant (for simplicity we consider a trivial global topol- 
ogy). We note that the smoothing gauge transformation 
is non-local, and so a careful study of potential problems 
associated with the infinite- volume limit is required. 

In the overlap approach, while the real part of the 
effective action is gauge invariant by construction, the 
imaginary part is not. Again, we expect that gauge- 
noninvariant counter-terms will be needed, starting at 
some finite loop order. Potentially severe problems with 
the overlap approach were pointed out in ref. |^^. Ac- 
cording to our judgement, subsequent works (including 
in particular ref. fail to address the issues raised in 
ref. psi . Numerical evidence for the lack of gauge invari- 
ance (in the non-abelian case) has been found in ref |3^] 
(see fig. 1 therein). 



V. THE REDUCED MODEL 

Returning to our approach, we consider in this section 
the limit of a vanishing gauge coupling. Since l/gg is 
the coefficient of the plaquette action, the go — limit 
constrains the lattice gauge field to the trivial orbit. 

The theory defined by (70 = limit is called the re- 
duced model. If we use the vector picture, the reduced 



in the 



model is obtained by substituting U. 
lattice action. The lattice gauge field measure Y[ J dUx,fi 
is replaced by J d(j)x ■ (Alternatively, starting from the 
Higg s pict ure that already involves both Ux,^ and (j)x (see 



Sect. [IB ), one obtains the reduced model by simply set- 
ting Ux,,, = /.) 

In the weak gauge-coupling limit, the transversal 
modes are perturbative at the lattice scale. Many impor- 
tant features, including the fermion spectrum, are deter- 
mined by the dynamics of the longitudinal modes. The 
utility of the reduced model is that it allows us to study 
the longitudinal dynamics in isolation, without making 
any a-priori assumption. The reduced model accounts 
for dynamical situations ranging from a divergent longi- 
tudinal correlation length, as in our approach, down to 
a very short correlation length. In this section we study 
a prototype reduced model. The entire dynamical range 
is realized in different regions of its phase diagram. In 
Sect. VI we discuss the effects of the longitudinal dynam- 
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ics on the fermion spectrum, first in general terms and 
then in our approach. 



A. The phase diagram 

In this subsection (and the next one) we discuss the 
reduced model related to t he la ttice action of Sect. 11 B , 
for (f)x £ U(l)- Since eq. (2^) is written in the Higgs 
picture, the reduced model is obtained by setting Ux,^, = 
1. This leads to the action 



5^ = E(-«' 



(5.1) 



Our first task is to derive the mean-field phase diagram 
in the (k, K)-plane. This phase diagram is in fact generic, 
and pertains also to the more relevant theories defined in 
Sect. and in //. 

Let us first consider the ordinary VEV, v, as an order 
parameter (or the staggered VEV vam)- An additional 
order parameter will be introduced shortly. By definition 



VAM 



(5.2) 
(5.3) 



where = (—1)^'' ''. For k = 0, we recover the fa- 
miliar non-linear sigma model. On the K-axis there is 
a symmetric (PM) phase for |k| < Kc, a ferromagnetic 
(FM) phase for k > Kc, and an antiferromagnetic (AM) 
phase for k < —Kc- The field redefinition (jj^ — > ex't'x 
maps K to —K, thus implying a symmetry of the K-axis. 

We now extend the discussion to the full (k, K)-plane. 
Mean-field approximation in c?-dimensions yields the fol- 
lowing equation for the FM-PM line 



K + {2d + l)k — Kc 



FM-PM. 



(5.4) 



The equation for the AM-PM line is 



K + {6d — 1)k 



AM-PM. 



The FM-PM and AM-PM transitions are continuous. 
The symmetry of the K-axis extends to k ^ 0. Under the 
field redefinition (j)^ — > Sxff'x, the point (k, k) is mapped 
in four dimensions to (k, — k — 32k). This implies that 
the linear equation 



16k = 0, 



(5.6) 



defines a symmetry line of the phase diagram. The FM- 
PM and AM-PM lines meet in the second quadrant, at 
the point (— -^j-!^^) on the symmetry line. It can be 
shown that, beyond this point, the symmetry line is a 
first-order transition line separating the FM and AM 
phases |3l| . 

In condensed matter physics, it is well-known that spin 
models with competing interactions tend to develop a 
ground state that breaks translation and rotation invari- 
ance. If a small antiferromagnetic interaction is added to 



a dominant ferromagnetic one, the spin orientation of the 
ground state will rotate slowly with a wave vector 
(see ref. 62] for a recent review). 



In the reduced model defined by eq. (5J^), competing 
interactions occur when k and k have opposite signs. In 
order to look for a similar phenomenon, we introduce the 
mean-field ansatz 



(5.7) 



We assume < < 27r. A non-zero signals the spon- 
taneous breaking of translation and rotation invariance. 
More precisely, translation invariance is broken in the 
direction defined by g^, but it remains unbroken in the 
transversal directions. In condensed matter, phases with 
a non-zero are known as helicoidal- ferromagnetic ones. 
Here, the helicoidal-ferromagnetic phase of the reduced 
model is the go — boundary of the FMD phase of the 
full theory (see Sect. |ll C| ), and the name FMD will be 
used both in the full theory and in the reduced model. 

A simple mean-field method is based on a factorized 
probability measure (see e.g. ref. ^^). For the FM-PM 
transition, one can use the following factorized probabil- 
ity measure 



Vo{0x) 



1 + 2ucos(6'j;) 
2^ 



(5., 



In order to accommodate a non-zero g^^, we generalize 
this to 



Vi9x 



1 + 2i;cos(6>3; - g^x^) 
2^ 



(5.9) 



One has (l)p — 1 and (e*^"=)p = ue'^^'^*' (in agreement 
with eq. (qJ)). Because of its factorized nature, the g^- 
dependence oiV{6x) affects only the internal energy, but 
not the entropy. Introducing the notation S'jj — Tlx , 
we find 



(5.5) (H)p = (1 - v^){Sn -f 72k) -f {nTiq^) + nT^{q^)) 



(5.10) 



where the function T is defined in eq. (2.8). While the 
choice of the factorized probability measures eqs. ( |5.8| ) 
and (p. 91) is somewhat arbitrary, the internal en- 
ergy (pTO) is a universal feature of any mean- field ap- 



proximation for S'^. 

If we consider a point in the phase diagram well to the 
right of the FM-PM line, the value of v is finite in lattice 
units. Making the self-consistent assumption that is 
small, the location of the FM-FMD transition can be de- 
termined by minimizing the internal energy with respect 
to g^. Remarkably, the g^j-dependent part of the internal 
energy coincides with the classical potential (2.7), if we 
make the identification goA^ ^ g^. Consequently, there 
is complete agreement between the mean-field properties 
of the FM-FMD transition in the reduced model, and the 
classical properties of the FMD transition in the 5o 7^ 
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theory. The mean-field location of the FM-FMD tran- 
sition is K = 0. For K > one is in the FM phase, 
whereas for k < one is in the FM D pha se. Close to 
the FM-FMD line, is given by eq. (|2.10[) where go^^ 
is replaced by (after this replacement go drops out, 
and one has = — k/2k for a small negative k). The 
FM-FMD hnc ends when it hits the FM-PM hne. The 
multi-critical point where the FM, FM and FMD phases 
meet is known as a Lifshitz point ||3^]. Its mean- field 
value is (^,0). (Lifshitz points exhibit rich critical be- 
haviour. This was discussed recently in a field theoretic 
context by J. Kuti 

The remaining features of the mean-field phase dia- 
gram are as follows (see FIG. 1). The transformation 
(f'x ixff'x maps q^ to q^ + tt. This implies the ex- 
istence of a second FMD region (and another Lifshitz 
point) below the symmetry line in the fourth quadrant. 
The PM-FMD line, separating the paramagnetic phase 
from the FMD phase, can be determined by first mini- 
mizing the (internal) energy with respect to q^, and then 
the (free) energy with respect to v. The FM phase oc- 
cupies a bounded region in the phase diagram. The two 
FMD regions above and below the symmetry line belong 
to a single FMD phase. The PM-FMD line hes on an 
ellipse, with the FM-PM and AM-PM lines tangent to it 
at the two Lifshitz points. A more detailed mean-field 



calculation will be presented elsewhere |31|. 

The mean-field ansatz ( ^.7| ) implies the simultaneous 
breaking of the internal U(l) symmetry, as well as of 
rotation and translation invariance. We believe that one 
cannot break translation invariance without at the same 
time breaking an internal symmetry. However, one can 
conceive of a phase (denoted PMD) where only rotation 
symmetry is broken, while the internal symmetries as 
well as translations are unbroken. The order parameter 
for a PMD phase of the reduced model is the expe ctati on 
value of the composite vector field (compare eq. (2.18)) 



1 



h.c. 



(5.11) 



At the moment, however, we have no evidence for a PMD 
phase. 

The phases of the reduced model are depicted in TA- 
BLE I. The order parameter vh is defined as = 
(0x6"^'^'''^^). Note that in the special case ~ 
(0,0,0,0) (q^ = (tt, TT, TT, tt)), Vh coincides with v (vam)- 
In this paper we usually do not distinguish between v 
and Vh, since the correct meaning can be understood 
from the context. However, this distinction is important 
in numerical simulation. The value of q^, to be used in 
the measurement of vh, can be determined for example 
by measuring (^(t)\.(t)x+fi) and extracting its phase. 

The relation between the {k, K)-phase diagram of the 
reduced model and the (it, k, (7o)-ptiase diagram of the 
full theory is the following. In the U(l) case, the sym- 
metric (PM) phase is the boundary of a Coulomb phase, 
and the broken (FM or AM) phase is the boundary of a 
Higgs phase. In the non-abelian case, the PM, FM and 



AM phases correspond to the boundary of a single Higgs- 
confinement phase. Finally, the helicoidal-ferromagnetic 
(FMD) phase of the reduced model is the boundary of 
the FMD phase of the full theory. 

In the go ~* limit of the full theory, we find (approx- 
imately) massless gauge bosons close to the FM-FMD 
line, as well as close to the PM-FM line, and in the entire 
PM phase. An interesting observation is that, even with- 
out fermions, if we want to study a Lorentz gauge-fixed 
Yang-Mills theory on the lattice, then the appropriate 
critical line is the FM-FMD line. The reason is that, in 
order to keep the longitudinal kinetic term in the tree- 
level action, k must scale like l/^g (see Sect. ||). In the 
large-K limit, a PM phase does not exist, and criticality 
can only be achieved by approaching the FM-FMD line. 



B. The weak-coupling expansion 
in the reduced model 

In view of the properties of the weak-coupling ex- 
pansion in the full theory (Sect. ||), and in particular 
eq. ( 2.12 ), one expects that 1/k will play the role of a 
coupling constant in the reduced model. We will now 
demonstrate this explicitly. We do not carry out here any 
detailed calculations that require the full lattice Feynman 
rules. Therefore, we work in the continuum approxima- 
tion, i.e. we extract from the lattice action the marginal 
and relevant terms, that control the critical behaviour in 
the vicinity of the gaussian critical point 1/k — Q. 

The weak-coupling expansion is facilitated by expand- 
ing around a broken symmetry vacuum. We first intro- 
duce the Goldstone boson (GB) field 0^ via (j)x = e'^=". 
(This classical expansion is consistent with the mean- 
field ansatz ( ^.7[ ), because f — > 1 for k ^ oo.) Rescaling 
9 {l/^/2k) 9 we find, in the continuum approximation, 
the following GB lagrangian 



^d^9d^9+^{d^9d^6 



(5.12) 



Eq. ( 5.12| ) is valid on the FM side of the transition line. 
On the FMD side, one first looks for the classical vacuum 
by assuming 9 — q^x^ and minimizing for q^. The result 
is the same as in the mean-field approximation. The 
weak-coupling expansion on the FMD side is then defined 
via 9 — > q^x^ + {l/^/2k) 9. 

We will consider here only the FM side. Taking the 
Fourier transform of the bilinear part of the lagrangian, 
we find the following GB propagator 



1 2 



(5.13) 



where mp = k/k. Analytical continuation to Minkowski 
space shows the existence of a positive-residue pole at 
p\i = 0, and a negative-residue (ghost) pole at p\j = tUq. 
These poles merge into a quartic singularity in the limit 
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0. The continuum limit of the reduced model 



is therefore not unitary. (In Sect. |VIB| we discuss the 
interaction of the GB field with fermions. The crucial 
requirement is that the non-unitary GB sector, which 
accounts for the two unphysical polarizations of the gauge 
bosons, will decouple from the fermions in the continuum 
limit. This decoupling is discussed in detail in ref. pGf.) 

Because of the quartic kinetic term in the GB la- 
grangian, the canonical dimension of the GB field 9{x) is 
zero. The GB lagrangian is invariant under the shift sym- 
metry 9{x) —^ 9{x) + const, which forbids the appearance 
of non-derivative terms under renormalization. In addi- 
tion, the GB lagrangian is inv arian t under the discrete 
symmetry 9{x) — > —9{x). Eq. (5.12) is the most general 
renormalizable lagrangian allowed by these symmetries. 

What marks the Feynman rules of the GB lagrangian, 
is that one derivative acts on every line attached to a 
vertex. The derivatives acting on the two ends of each 
internal line effectively cancel one factor of 1/p^ in the 
propagator. The result is that the UV power counting of 
the GB model is the same as in an ordinary theory. 
If we ignore the vector index carried by the partial deriva- 
tives, one can match each term in the GB lagrangian with 
a corresponding term in the A^"* lagrangian according to 
the rule 96* ^ <&. 

In the A*!''* theory, at the one loop level only the mass 
term is renormalized, but not the kinetic term. By anal- 
ogy, in the GB lagrangian only (9^6*)^, but not (□6')^, is 
renormalized at the one-loop level. The induced one-loop 
(i9^6')^-term has a positive sign. This implies 



Kc.l{k) = -c+0 (k ^) 



(5.14) 



(Note that the coefficient of {d^OY in eq. ( 5.12 ) is k/k.) 
The dimension of the positive constant c is two. Its nu- 
merical value, which is 0(l/a^), has to be determined 
by a lattice calculation. Finally, as in the A$'' theory, 
the one-loop beta-function is determined by the vertex 
renormalization, and is found to be positive. Explicitly 



IGtt^ 



(5.15) 



where A is the UV cutoff (the inverse lattice spacing in 
a lattice calculation). 



The IR divergences of massless Goldstone bosons lead 
to the restoration of continuous symmetries in two di- 
mensions Only symmetric observables, which are 
invariant under all the continuous symmetries, can have 
a non-zero value. There are theorems [p8|-^ that guar- 
antee the IR-finiteness of the s ymme tric observables. 

Here the quartic propagator ( ^.16 ) does not character- 
ize a whole phase, but only the FM-FMD line itself. The 
order parameter vh (or v) dips close to the FM-FMD 
line, and vanishes on that line in several (may be in all) 
interesting cases The theorems on the finiteness 

of symmetric observables, in particular ref. gener- 
alize to four dimensions. Also, as in two dimensions, 
the predictions of the weak-c oup ling expansion are often 
valid, if interp reted carefully |^1|]. This will be important 



in Sect. VI B 



D. Realistic reduced models 

The key features of the simplified reduced model stud- 
ied in this section extend to the realistic reduced mod- 
els, d efine d from the (gauge-fixing and ghost) action of 
Sect. [IE for the non-linear gauge, or the action of // for 
the linear gauge. This includes the qualitative structure 
of the phase diagram and in particular the FM and FMD 
phases, the gaussian critical point at k = oo, and the IR 
divergences on the FM-FMD line. 

Particularly interesting are the critical FM-FMD the- 
ories in the reduced models that correspond to the linear 
gauge d ■ A = Q. In the abelian case, the linear-gauge re- 
duced model leads to a free theory with a l/(p^)^ propa- 
gator. The properties of this critical theory are analogous 
to the spin- wave phase of a two-dimensional abelian the- 
ory. This will be discussed in detail elsewhere pl| , p6[ . 
The critical theory for a non-abelian gauge group was in- 
vestigated by Hata in the continuum approximation. 
His main result is that, like non-abelian sigma models in 
two dimensions, these four-dimensional non-linear mod- 
els are asymptotically-free. It will be interesting to inves- 
tigate the significance of this result for the construction of 
gauge-fixed non-abelian lattice theories via our approach. 



VI. FERMIONS IN THE REDUCED MODEL 



C. Infra- Red divergences of the critical theory 

If we tune k to Kc./.(k), the quadratic kinetic term in 
eq. ( 5.13| ) vanishes, and the renormalized GB propagator 
reads 



(5.16) 



where Z accounts for the wave-fuction renormalization. 
This quartic propagator leads to IR divergences in four 
dimensions, like massless bosons with an ordinary kinetic 
term do in two dimensions. 



In a manifestly gauge invariant theory like QCD, the 
fermion spectrum can be read off from the lattice ac- 
tion by going to the free field limit = 0. Here, the 
fermion action is not gauge invariant (in the vector pic- 
ture), and the limit go = gives rise to an interact- 
ing theory, namely, to the reduced model. We identify 
the elementary fermions of a general lattice gauge theory 
with the independent fermionic massless poles of the as- 
sociated reduced model. (If the fermion action is gauge 
invariant, any (/)2;-dependence of its reduced-model form 
can be eliminated by a field redefinition.) It is justified to 
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determine the matter spectrum by setting go = 0, since, 
in a scahng region, the transversal degrees of freedom are 
perturbative at the lattice scale. 



A. The robustness of the No-Go theorems 

Let the gauge field belong to a Lie group G. By con- 
struction, the associated reduced model has a global G- 
symmetry, denoted Gl, that acts on (px by left multipli- 
cation. (The reduced model is obtained from the vector 



picture via C/^, 



and the product 



invariant under left multiplication. Notice also that the 
gauge-invariant Higgs picture can be obtained by gaug- 
ing the Gl symmetry of the reduced model.) Now, we 
demand the existence of massless vector bosons in the 
scaling region, which can be identified with the gauge 
bosons of the target continuum theory. These vector 
bosons couple to the Noether current associated with the 
Gl symmetry. Thus, assigning the fermions to represen- 
tations of Gl determines whether the continuum limit is 
chiral or vector-like. 

In previous chiral fermion proposals, it was usually 
attempted to take the continuum limit in a symmetric 
phase, where the Gl symmetry is not broken sponta- 
neously. (Physical gauge invariance is restored dynam- 
ically in a symmetric phase, when we consider the full 
50 7^ theory. This means that there are no light un- 
physical states, whose decoupling in the continuum limit 
requires fine-tuning. Since the VEV of the (px field is 
zero, the physics in a symmetric phase is more easily ac- 
counted for in the gauge- invariant Higgs picture.) In a 
symmetric phase, the fluctuations of the <j)x field are usu- 
ally not controlled by any small parameter. As a result, 
non-perturbative methods had to be invoked in order to 
determine the fermion spectrum. Where available, it was 
always found that the true fermion spectrum is vector- 
like (see ref. @,|,|| for details). 

We have discussed this phenomenon in ref. and 
argued that it has a simple physical explanation. Here 
we can only outline the key considerations leading to this 
conclusion, and we refer the reader to ref. [^j^ for the de- 
tails. One starts with the observation that, in a symmet- 
ric phase of the reduced model, there are generically no 
massless scalars. Therefore, the only massless particles 
(if any) are fermions. (It could be that no massless 
fermions are present unless a mass term if fine-tuned. 
Since we are in symmetric phase, a massless fermion 
obtained by fine-tuning is necessarily a Dirac fermion.) 
Now, in four dimensions, there are no renormalizable in- 
teractions involving only fermion fields. The continuum 
limit defined by a generic point inside a symmetric phase 
is therefore a theory of free massless fermions (if it is 
not empty). One can then construct an effective lat- 
tice hamiltonian for the fermions, that satisfies all the 
assumptions of the Nielsen-Ninomiya theorem. (The ef- 
fective hamiltonian is defined as the po = limit of the 



inverse of a suitable two-point function.) We refer here in 
particular to the analytic structure near the zeros of the 
effective hamiltonian, and to the existence of a smooth 
interpolation throughout the rest of the Brillouin zone. 
This leads to the conclusion that the fermion spectrum is 
vector-like in a symmetric phase, provided the underly- 
ing theory is local. (In the case of a non-local theory one 
expects violations of unitarity and/or Lorentz invariance, 
see ref. for references to the original literature.) 

This impasse extends, by continuity, to the fermion 
spectrum on any phase transition line that separates a 
symmetric phase from a broken phase. In particular, 
even though the gauge boson mass vanishes on the PM- 
FM line, we do not expect to find a chiral gauge theory 
by taking the continuum limit at the PM-FM line. The 
fermion spectrum will be vector-like if the PM-FM line 
is approached from the PM phase. If we approach the 
PM-FM line from the FM phase, we can only obtain a 
mirror fermion model p^ , but we cannot decouple the 
unwanted mirror fermions. 



B. Evading the No-Go theorems 

Let us now investigate what changes when the contin- 
uum limit is taken at the FM-FMD line. We will consider 
the simplest case, namely a U(l) gauge group with the 
gauge-fixing action pertaining to the linear gauge, cf. II. 
As mentioned in Sect. VD| , the properties of the critical 
FM-FMD theory are analogous to the spin- wave phase of 
a two-dimensional abelian theory. 

We go from (the vector picture of) the full theory to 
the reduced model according to the rule Ux.^i. 4)\.(t)x+ji- 
The fermion action eq. (pj) becomes (we use the two- 
component notation) 



S'f = E(x0^)^mA;.(</>x) - J (xux + fi.' 



(6.1) 



The fermion variables in eq. (6.1) are neutral wit)! respect 
to Gl- If, instead, we use the c/iarged variables Xc — 'PXi 
the fermion action reads 

= ^xcCT^A^Xc - j(('/'^xc)n(0^xc) + fi.c. 



(6.2) 

According to the rules of the weak-coupling expansion 
(see Sect. VB), the tree-level fermion action is obtained 
by substituting the classical vacuum (jjx — 1. Using 
eq. (^) for definiteness, we get 



.XcDXc 



(6.3) 



In the limit w ~ 0, only the kinetic term Xc A^x^ 
is left. Thus, the w = action exhibits the infamous 
doubling, with sixteen Weyl fermions altogether. (Each 
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fermion is associated with a point in the Brillouin zone, 
whose lattice momentum components are equal to either 
or TT.) Since we take w — 0(1), the MW term elimi- 
nates the doublers, and the pole in the tree-level fermion 
propagator describes a single Weyl field (see Sect. ill). 
Had we started from the fermion action written in 



terms of the neutral variables (eq. (6.1)), the substitu- 
tion (hj. — 1 would lead to a tree-level action identical 



to eq. (S.3), but with the neutral field x replacing the 
charged field Xc- Now, deep in the FM phase this makes 
no difference, because the Gl symmetry is broken any- 
way by the 0a;-VEV, which is 0(1) in lattice units. How- 
ever, the Gl symmetry is restored right on the FM-FMD 
line It is therefore a meaningful (and important) 

question to ask what are the GL-quantum numbers of 
the massless fermions. 

The fact that one cannot simply read off the quan- 
tum numbers of one-fermion states from the tree-level 
action, is a consequence of the IR-divergent nature of 
the GB propagator, l/(p^)^. The way to proceed is 
to examine a family of fermionic two-point functions 
Tn = (('/'"x) (X'/'^"))- Assuming all mass parameters 
have been tuned to their critical values, r„ will in gen- 
eral contain terms proportional to {'f)~^ \o^{p^) for any 
k. The presence of logarithmic terms (which typically 
lead to power law corrections when summed over all or- 
ders) means that the operator (jf^x does not create a one- 
particle state 1 41 4^ . Only when there are no logarithmic 
corrections do we have a simple massless pole, and the 
quantum numbers of the intermediate one-fermion state 
must coincide with the quantum numbers of the interpo- 
lating fermion field. 

The fermion spectrum in the reduced model can 
be studied in detail using the weak-coupling expan- 
sion. While the actual calculations require a substan- 
tial amount of work, the conclusions are robust, as they 
really depend on universal properties of the low energy 
effective (continuum) lagrangian. A one-loop calculation, 
which is also supported by numerical simulations, will be 
presented elsewhere Here we will list the key results, 
as they apply to the MW fermion action. 

• Logarithmic terms are absent only for n = 1, 
namely in the two-point function (xcXc)- Con- 
sequently, the massless fermion has the quantum 
numbers of the Xc field. The latter is charged, which 
means that the Xc Weyl fermion will couple to the 
transversal gauge field, when the latter is turned 
on. 

• As can be expected on general grounds (see 
Sect. [H), divergent Majorana-like mass terms are 
induced at the one-loop level. These must be can- 
celled by suitable counter-terms, to maintain the 
masslessness of each chiral fermion. 

• When the Majorana-like mass counter-terms are 
tuned to their critical values, the unphysical GB 



field decouples from the Xc-fermions. The contin- 
uum limit is a direct product of (in general) several 
free theories, one associated with the unphysical 
GB field, and one associated with every species of 
chiral Xc-fermions. 

This establishes an agreement between the predictions of 
the weak-coupling expansion in the full theory and in the 
reduced model, thus supporting the consistency of our 
approach. The properties of the reduced model are true 
for an arbitrary fermion spectrum, and this is consistent 
with the vanishing of the anomaly in the absence of a 
transversal gauge field. 

In com pariso n with previous chiral fermion proposals 
(see Sect. VIA) we note two key differences that allow us 
to escape from a similar impasse. First, one may worry 
that the need to tune mass counter-terms may indicate 
that we got the wrong spectrum {e.g. Dirac instead of 
Weyl fermions). Now, when the (Majorana-like) mass 
counter-terms are not tuned to their critical values, the 
fermions remain coupled to the IR-singular GB field in 
the low-energy limit. Due to potential IR divergences, it 
is not at all clear that (massive) one-fermion states could 
be consistently defined in this case, nor that such states 
would have well-defined GL-quantum numbers. The off- 
critical theory remains to be investigated in the future. 
However, in view of the above IR subtleties, the general 
conclusion is that by considering the role of fermion mass 
perturbations, one does not end up with an argument 
against the existence of a chiral spectrum at the critical 
point. (See also the discussion of fermion mass counter- 
terms in Sect. [II.) 

The other key difference is that the continuum limit is 
now taken at the phase transition separating two broken 
phases of the reduced model. Off the FM-FMD line (on 
both sides) the Gl symmetry is broken spontaneously, 
and all asymptotic states do not have well-defined Gl- 
quantum numbers. On the FM-FMD line itself, the Gl 
symmetry is restored, and the question arises whether 
we do not run into the same old confiict with the No- 
Go theorems. The answer is contained in the analytic 
structure discussed above. Thanks to the presence of 
the highly IR-singular GB field, a zero in the inverse 
propagator does not necessarily imply the existence of a 
one-fermion state with the same quantum numbers. As 
an illustration, consider the four-component unit-charge 
field ^^l)c^ whose left-handed component is PliPc — Xcj 
and whose right-handed component is PriPc = 4''^Xc- If 
we consider the inverse two-point function of we may 
erroneously conclude that it interpolates a massless Dirac 
fermion. In reality, only the left-handed channel of this 
inverse propagator has a simple zero ~ ^, implying the 
existence of a unit-charge left-handed fermion. In the 
right-handed channel, one finds a j^log(p^) correction in 
the one-loop approximation, which implies the absence 
of a right-handed fermion with the same charge. 
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VII. OPEN QUESTIONS 



In Sect. [IE, the criterion for fixing the counter-terms 



was to enforce BRST invariance (and, hence, unitarity) 
order by order. This perturbative prescription is incom- 
plete. Ultimately, the counter-terms should be deter- 
mined by a non-perturbative method. To rigorously de- 
fine the continuum limit, one has to specify a trajectory 
in the Higgs (or Higgs-confinement) phase, that ends at 
the gaussian p oint qq = 1/k = on the FMD bound- 
ary. (See Sect. V A for the phase diagram.) In addition, 
one has to construct a BRST operator, and prove its 
nilpotency in the continuum limit. Enforcing BRST in- 
variance should also lead to the restoration of full SO (4) 
invariance, because the marginal S0(4)-breaking opera- 
tors violate the BRST symmetry too. 

We comment that similar problems are encountered 
in lattice QCD with Wilson fermions, where the axial- 
flavour symmetries are broken on the lattice, in anal- 
ogy with the BRST symmetry in our gauge-fixing ap- 
proach. When using Wilson fermions, tuning is required 
not only at the level of the lattice action, but also in 
the construction of renormalized operators with well de- 
fined axial-flavour transformation properties [^5| . This 
is analogous to the problem of deflning BRST-invariant 
operators in our gauge-flxing approach. (In QCD, the 
flne-tuning problem can be solved using domain-wall 
fermions Q-Q. Whether a similar solution exists for 
the tuning problem in the gauge-flxing approach, is an 
interesting question.) 

Our gauge-fixing formulation can be tested by applying 
it to asymptotically-free gauge theories which are not chi- 
ral. In particular, in the absence of fermions, one should 
study whether the confining behaviour and the mass gap 
of Yang-Mills theories are reproduced. One possibility is 
that the FMD transition becomes weakly first-order due 
to non-perturbative effects. This scenario is favourable, 
at least from the point of view of numerical simulations. 
Another possibility is that the correlation length of the 
vector field strictly diverges at the FMD transition, al- 
ready for go ^ 0. In this case, a consistent continuum 
limit may exist provided all the massless excitations are 
unphysical. 



VIII. CONCLUSIONS 

In a regularized chiral gauge theory, the longitudi- 
nal modes of the gauge field couple to the fermions. 
Before the regularization is removed, there are viola- 
tions of gauge invariance even if the fermion spectrum 
is anomaly-free. When we use the lattice regularization, 
the longitudinal modes should decouple in the contin- 
uum limit, but it may be too much to expect for exact 
decoupling when the lattice spacing is still finite. 

The gauge-fixing approach aims to decouple the lon- 
gitudinal modes in the continuum limit. In this paper 



we have discussed how the gauge-fixing approach may 
be realized, thus making the first step of a systematic 
investigation of the gauge fixing approach. We have con- 
structed a lattice gauge-fixing action that has a unique 
classical vacuum. The gauge-fixing action contains a 
longitudinal kinetic term, and leads to a renormalizable 
weak-coupling expansion, which is valid even if the lattice 
fermion action is not gauge invariant. We have argued 
that the continuum fields, needed to describe the scal- 
ing behaviour, are in one-to-one correspondence with the 
poles of the tree-level lattice propagators. This should 
accommodate any consistent theory, including anomaly- 
free chiral gauge theories. 
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FIG. 1. Mean-field phase diagram. See TABLE I for the 
definition of the various phases. 



TABLE L Phases of the reduced model. The entries indi- 
cate which order parameters are non-zero in each phase. 



phase 


V 


VAM 


VH 




PM 


no 


no 


no 


no 


FM 


yes 


no 


yes 


no 


AM 


no 


yes 


yes 


no 


FMD 


no 


no 


yes 


yes 


PMD(?) 


no 


no 


no 


yes 



16 



